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Comments on 
"Modified negative binomial description of the 
multiplicity distributions in lepton-nucleon scattering" 

S. HegyiQ 

Research Institute for Particle and Nuclear Physics, 
H-1525 Budapest 114, P.O. Box 49. Hungary 

In two recent papers the present author introduced a generalization of the 
Negative Binomial Distribution [1,2]. It was obtained by extending the va- 
lidity of the asymptotic KNO scaling form of the NBD to the /x > powers 
of the scaling variable x = n/n. With this modification the asymptotic KNO 
function becomes 



A^^x^^'-^exp(-[Ax] 



(1) 



which is the generalized gamma distribution. In pQCD /i = (1 — 7)^^ with 
7 oc ^/al being the QCD multiplicity anomalous dimension. Poisson trans- 
forming f{x) one obtains the cited generalization of the NBD whose analytic 
form is given by 
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for /i > 1 
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" n\T{k) 

where 6 = nr{k)r~^{k + l/fi) and H(-) is the Fox generalized hypergeometric 
function or H-function. For /i = 1 the NBD is recovered. For /i > 1 P„ is 
not infinitely divisible and the factorial cumulants of the distribution exhibit 
nontrivial sign-changing oscillations. Eqs. (2)-(3) involve as special and lim- 
iting cases the Poisson transform of numerous classical probability densities 
such as the Chi, WeibuU, Maxwell, Pareto and Lognormal to mention but a 
few. Since the name Generalized NBD is already in use (not rarely to denote 
different discrete distributions) we shall call the above H-function extension 
of the Negative Binomial as HNBD for short. 
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In ref. [2] the HNBD was used to analyse the charged particle multi- 
plicity distributions measured by the HI Collaboration in deep inelastic e+p 
scattering at HERA [3]. Similar analysis was performed in ref. [4] using an- 
other 3-parameter generalization of the NBD, to so-called Modified NBD. 
The generating function G{z) — Z^^q -fn^" of the MNBD is given by 

= with r^A + n/k (4) 

which yields for A = the generating function of the NBD. The analytic 
form of the MNBD is usually written as 

During the past few years the MNBD was successfully used to describe the 
multiplicity distributions measured in different collision processes. For the HI 
data the quality of the MNBD fits reported in [4] are again very satisfactory. 
The author of ref. [4] made a comparison of the MNBD and HNBD results. 
It was found that the x^/d.o.f. values obtained by the one-parameter MNBD 
fits are significantly smaller than those obtained by the two-parameter HNBD 
fits. This seems to indicate that the HNBD is less satisfactory in reproducing 
the HI data than the MNBD. But the reason behind the better quality of 
MNBD fits is different. 

In ref. [2] the HNBD was fitted to the corrected HI data with variable 
parameters n and The shape parameter was fixed at A; = 1 (WeibuU case). 
As is clearly stated in [2] the usually less than 5% relative errors of the best-fit 
parameters are only statistical, the systematic uncertainties of the HI data 
were not taken into account in the analysis. The MNBD was examined in 
ref. [4] with fixed shape parameter k = 7 and variable A. The value of n was 
taken from ref. [3]. Although it is not stated explicitly, the 15-25% relative 
errors of the single fit parameter A indicate that the MNBD analysis was 
carried out using both the statistical and systematic uncertainties of the HI 
data. The comparison of the quality of fits obtained in refs. [2] and [4] is 
therefore misleading. 



2 



To clarify the situation the present author repeated the HNBD analysis 
with the inclusion of both the statistical and systematic errors of the HI 
data. Following ref. [4] n was fixed at its experimental value and only ix 
was taken as free parameter (A; = 1 as earlier). The comparison of the 
MNBD and HNBD fits is shown in Table 1. We can see that the HNBD 
easily reproduces the quality of MNBD fits. Of course it is not particularly 
meaningful to speculate which distribution is better; neither of them can be 
ruled out by a x^-analysis of the HI data tabulated in [3]. As a more stringent 
test, it would be of interest to compare MNBD and HNBD fits with the use 
of the full covariance matrix which makes possible to take into account the 
correlations between adjacent multiplicities. 



W (GeV) 


-A 


xV d.o.f. 




XV d.o.f. 


80 -f- 115 


0.216 ±0.034 


2.7/18 


4.811 ±0.392 


0.9/18 


115^ 150 


0.233 ±0.036 


10.6/21 


4.238 ±0.277 


2.2/21 


150^ 185 


0.187 ±0.049 


4.6/22 


4.011 ±0.288 


1.9/22 


185 ^ 220 


0.264 ±0.056 


11.3/23 


4.021 ±0.334 


1.5/23 



Table 1. Comparison of the MNBD fits, taken from ref. [4], and the HNBD 
fits for the HI multiplicity data corresponding to the pseudorapidity interval 
1 < ?7* < 5. The errors quoted are statistical and systematic. 
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